home *** CD-ROM | disk | FTP | other *** search
/ Info-Mac 4 / Info_Mac IV CD-ROM (Pacific HiTech Inc.)(August 1994).iso / Science / MAL ©P.f.Howden 1⁄1⁄89 / LOGLOGFIT < prev    next >
Text File  |  1987-07-27  |  1KB  |  19 lines
  1. 5 CLS:PRINT"LOG-LOG DATA FIT":CLEAR:DEFDBL A-Z:INPUT "MENU:0=STRAIGHT LINE: 1=CURVE:=";E
  2. 7 IF E=1 THEN PRINT"THE FUNCTION WILL BE: LN Y=A*(LN X)^N+B*LN X+C":GOTO 10
  3. 8 PRINT"THE FUNCTION WILL BE: LN Y=A*LN X+B"
  4. 10 PRINT"INPUT DATA POINTS:":INPUT "X1=";X1:INPUT "Y1=";Y1:INPUT "X2=";X2:INPUT "Y2=";Y2:IF E=1 THEN 50
  5. 15 A=(LOG(Y1)-LOG(Y2))/(LOG(X1)-LOG(X2)):PRINT"FUNCTION IS: LN Y=";A;"*LN X+(";
  6. 20 B=LOG(Y1)-A*LOG(X1):PRINT B;")":PRINT
  7. 25 INPUT "TO CHECK,INPUT ANY DATA POINT, X=";X:PRINT"CORRESPONDING Y=";EXP(B+A*LOG(X)):GOTO 25
  8. 50 INPUT "X3=";X3:INPUT "Y3=";Y3
  9. 55 INPUT "DEGREE N (TRY AN INTEGER >=2)=";N
  10. 60 A=(LOG(Y1)-LOG(Y2))*(LOG(X2)-LOG(X3))-(LOG(Y2)-LOG(Y3))*(LOG(X1)-LOG(X2))
  11. 65 A=A/(((LOG(X1))^N-(LOG(X2))^N)*(LOG(X2)-LOG(X3))-((LOG(X2))^N-(LOG(X3))^N)*(LOG(X1)-LOG(X2)))
  12. 70 B=(LOG(Y1)-LOG(Y2)-A*((LOG(X1))^N-(LOG(X2))^N))/(LOG(X1)-LOG(X2))
  13. 75 C=LOG(Y1)-B*LOG(X1)-A*(LOG(X1))^N
  14. 80 PRINT"FUNCTION IS:  LN Y=";A;"*(LN X)^";N:PRINT"+(";B;")*LN X+(";C;")"
  15. 85 PRINT:INPUT "TO CHECK,INPUT ANY DATA POINT, X=";X:PRINT"CORRESPONDING Y=";EXP(A*LOG(X)^N+B*LOG(X)+C):INPUT "CHANGE DEGREE (Y/N) ,else continue checking";A$
  16. 90 IF A$="Y" THEN 55
  17. 95 GOTO 85
  18.          MACINTOSH LISTING
  19.